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Abstract 

This presentation starts with a brief review of our current picture of QCD phases, derived from 
lattice QCD thermodynamics and from models based on the symmetries and symmetry breaking 
patterns of QCD. Typical approaches widely used in this context are the PNJL and chiral quark- 
meson models. It is pointed out, however, that the modeling of the phase diagram in terms of 
quarks as quasiparticles misses important and well known nuclear physics constraints. In the 
hadronic phase of QCD governed by confinement and spontaneously broken chiral symmetry, in- 
medium chiral effective field theory is the appropriate framework, with pions and nucleons as active 
degrees of freedom. Nuclear chiral thermodynamics is outlined and the liquid-gas phase transition 
is described. The density and temperature dependence of the chiral condensate is deduced. As a 
consequence of two- and three-body correlations in the nuclear medium, no tendency towards a 
first-order chiral phase transition is found at least up to twice the baryon density of normal nuclear 
matter and up to temperatures of about 100 MeV. Isospin- asymmetric nuclear matter and neutron 
matter are also discussed. An outlook is given on new tightened constraints for the equation- 
of-state of cold and highly compressed matter as implied by a recently observed two- solar- mass 
neutron star. 

Introduction 

QCD phase diagram: visions and facts 

A theorist's sketch of the QCD phase diagram, plotted in the plane of temperature T and baryon 
chemical potential usually looks like the one shown in FigJTJ The hadronic phase, with quarks and 
gluons confined in mesons and baryons, is thought to be separated from the deconfined quark-gluon 
phase by a crossover transition at low baryon chemical potential \ib and at temperatures T of order 
Aqcd ~ 0.2 GeV. At a critical point this crossover supposedly ends and turns, at larger values of the 
baryon chemical potential, into a first-order phase transition. When the baryon chemical potential is 
further increased, various low-temperature color-superconducting phases are expected to emerge. 

At zero chemical potential, lattice simulations of full QCD thermodynamics with three quark flavors 
[H El El Hj do indeed suggest a crossover scenario for the temperature dependence of the chiral (quark) 
condensate, an order parameter spontaneously broken chiral symmetry, with a transition temperature 
T c around 170 MeV. Lattice QCD extrapolations to non-zero (real) chemical potentials are notoriously 
difficult, however. Extensions of the phase diagram to baryonic matter and various forms of quark 
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Figure 1: Schematic picture of the QCD phase diagram in terms of temperature T and baryon chemical 
potential /ig, sketching the deconfined quark-gluon phase at high temperature, the low-temperature and 
low-density hadronic phase with spontaneously broken chiral symmetry, and the domain of high-density 
(color-superconducting - CSC) phases featuring various forms of quark Cooper pairing. The dashed line 
illustrates a chiral crossover transition ending at a hypothetical critical point where a first-order phase 
transition boundary (solid line) begins. The first-order liquid-gas phase transition of nuclear matter is 
also indicated. 



matter rely so far entirely on models. Only at asymptotically large quark chemical potentials the 
picture becomes simple again. At Fermi momenta of the order of several GeV, the quarks at the Fermi 
surface can be treated using perturbative QCD. The gluon exchange quark-quark interaction is weakly 
attractive in color-antitriplett pairs. Such a weakly attractive interaction is sufficient to form Cooper- 
pair-like diquarks in various combinations so that color superconducting phases emerge [5]. However, 
the more quantitative density scale at which such transitions take place cannot be predicted. 

Basically all models of the QCD phase diagram (briefly summarized in Section 2) work with dynam- 
ical quarks as quasiparticles and do not take into account baryons and their correlations. An apparently 
generic result of such models is the fact that the first order transition line at low temperature meets 
the chemical potential axis at a quark chemical potential around ji q ~ 0.3 GeV, the typical constituent 
quark mass scale, corresponding to a baryon chemical potential \±b close to the nucleon mass. This, 
however, is the terrain of nuclear physics with its enormously rich and well explored phenomenology. 
No sign of a first-order chiral or deconfinement phase transition is visible in that area right in the center 
of the QCD phase diagram. The only known phase transition in nuclear matter is the one from a 
Fermi liquid to an interacting Fermi gas. This motivates the present essay as a proposal for taking the 
well-known nuclear physics constraints more seriously into account (see Section 3). This presentation 
is in parts an update of a previous review [6]. It closes with a brief discussion of new constraints from 
neutron star observations on the equation of state at high density and low temperature. 
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2 Modeling the QCD phase diagram 



Confinement and spontaneous chiral symmetry breaking are governed by two basic symmetry principles 
of QCD: 

• The symmetry associated with the center Z{3) of the local SU(3) C color gauge group is exact 
in the limit of pure gauge QCD, realized for infinitely heavy quarks. In the high-temperature, 
deconfinement phase of QCD this Z{3) symmetry is spontaneously broken, with the Polyakov 
loop acting as the order parameter. 

• Chiral SU(Nf) R x SU(Nf) L symmetry is an exact global symmetry of QCD with Nf massless 
quark flavors. In the low-temperature (hadronic) phase this symmetry is spontaneously broken 
down to the flavor group SU{Nf) v (the isospin group for Nf = 2 and the "eightfold way" for 
Nf = 3). As a consequence there exist Nf — 1 pseudoscalar Nambu-Goldstone bosons and the 
QCD vacuum hosts a strong quark condensate. 

These symmetries and symmetry breaking patterns serve as guiding principles for constructing 
models of the QCD phases. The basic question is about the intertwining of spontaneous chiral symmetry 
breaking with confinement. 

2.1 Chiral condensate and Polyakov loop 

The order parameter of spontaneously broken chiral symmetry is the quark condensate, (qq). The 
disappearence of this condensate, by its melting above a characteristic transition temperature scale, 
signals the restoration of chiral symmetry in its unbroken Wigner-Weyl realization. The transition from 
confinement to deconfinement in QCD is likewise controlled by an order parameter, the Polyakov loop. 
A non- vanishing Polyakov loop $ reflects the spontaneously broken Z(3) symmetry characteristic of the 
deconfinement phase. The Polyakov loop vanishes in the low-temperature, confinement sector of QCD. 

Two limiting cases are of interest in this context. In the pure gauge limit of QCD, corresponding 
to infinitely heavy quarks, the deconfinement transition is established as a first order phase transition 
with a critical temperature of about 270 MeV. In the limit of massless u and d quarks, on the other 
hand, the isolated chiral transition appears as a second order phase transition at a significantly lower 
critical temperature. This statement is based on calculations using Nambu - Jona-Lasinio (NJL) type 
models M> E] which incorporate the correct spontaneous chiral symmetry breaking mechanism but 
ignore confinement. The step from first or second order phase transitions to crossovers is understood 
as a consequence of explicit symmetry breaking. The Z{3) symmetry is explicitly broken by the mere 
presence of quarks with non-infinite masses. Chiral symmetry is explicitly broken by non-zero quark 
masses. It is then a challenging question whether and how the chiral and deconfinement transitions 
get dynamically entangled in just such a way that they finally occur within overlapping transition 
temperature intervals. 

Confinement implies spontaneous chiral symmetry breaking, but the reverse is not necessarily true. 
There is no a priori reason why the chiral and deconfinement transitions should be intimitely connected. 
Nonetheless this appears to be the case in lattice QCD computations with almost physical quark masses. 
Results from lattice QCD thermodynamics [lj with 2+1 flavors (at zero baryon chemical potential) 
gave a chiral crossover transition temperature of about 190 MeV. Recent improved computations j2j [3] 
indicate a shift of the chiral transition to somewhat lower temperatures, around 160 MeV, consistent 
with earlier lattice simulations [4j. Lattice QCD results for the temperature dependence of the Polyakov 
loop have remained remarkably stable over the years. Pure gauge QCD on the lattice, with full gluon 
dynamics but without quarks (or equivalently, with infinitely heavy quarks) has established a first-order 
deconfinement transition at a critical temperature T c ~ 270 MeV. Once light quarks are added, the Z{3) 
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Figure 2: Pure gauge and full QCD lattice results for the Polyakov loop as function of temperature pQ. 
The solid line is calculated in an advanced (non-local) PNJL model [15], to be discussed later in the 
text. 



center symmetry is explicitly broken by the presence of these quarks and the first-order deconfinement 
transition turns into a crossover, with a significantly reduced transition temperature around 200 MeV. 
This is shown in Fig j2] where Polyakov loop results from pure gauge lattice computations are compared 
to full QCD simulations with Nf = 2 + 1 quark flavors. 

2.2 The non-local PNJL model 

Insights concerning the issue of a possible intertwining of chiral and deconfinement transitions can 
be gained from a model based on a minimal synthesis of the NJL-type spontaneous chiral symmetry 
breaking mechanism and confinement implemented through Polyakov loop dynamics. This PNJL model 
fTUl [TT] is specified by the following action: 



r/3=l/T r r w 

S= dr d 3 4^-%^i) - -W($,T) . (l) 

Jo Jv 1 J 1 

It introduces the Polyakov loop, 

$ = N- 1 Trexp(i</>/T) , (2) 

with a homogeneous temporal gauge field, = 3 A 3 + 8 A 8 G £77(3), coupled to the quarks. The 
dynamics of $ is controlled by a Z(3) symmetric effective potential designed such that it reproduces 
the equation of state of pure gauge lattice QCD with its first order phase transition at a critical 
temperature of 270 MeV. The field (j) acts as a potential on the quarks represented by the flavor doublet 
(for Nf = 2) or triplet (for Nf = 3) fermion field ip. The Hamiltonian density in the quark sector is 

H = -i^(d • V + 74 m - ^ + V(^, ^) , (3) 

with the quark mass matrix rh and a chiral SU{Nf)i, x SU{Nf)n symmetric interaction V. 

Earlier two-flavor versions of the PNJL model [T0l[TTl[T2] have still used a local four-point interaction 
of the classic NJL type, requiring a momentum space cutoff to regularize loops. A more recent version 
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J, using a non-local interaction between quarks, does not need any longer an artificial NJL cutoff. It 
generates instead a momentum dependent dynamical quark mass, M(p), along with the non- vanishing 
quark condensate. A further extension to three quark flavors [14] includes a U(1)a breaking term 
implementing the axial anomaly of QCD. This term is constructed as a non-local generalization of the 
Kobayashi-Maskawa-'tHooft 3x3 determinant interaction. 

A basic relation derived in this non-local PNJL model is the gap equation generating the momentum 
dependent quark mass self-consistently together with the chiral condensate (^). Its form (written 
here for the two-flavor case) is reminiscent of the corresponding equation emerging in Dyson-Schwinger 
approaches to QCD: 

where m = m u = is the current quark mass, G is a coupling strength of dimension (length) 2 and 
C(q) is a momentum space distribution, the Fourier transform of which represents the range over which 
the non-locality of the effective interaction between quarks extends in (Euclidean) 4-dimensional space- 
time. This non-local approach permits to establish contacts with the high-momentum limit of QCD 
with its well-known behaviour M(p) oc —a s (p)(ipilj)/p 2 at p — >> oc. At p < 1 GeV the distribution C(p) 
is designed to follow lattice QCD and Dyson-Schwinger results, or it may be deduced from an instanton 
liquid model (see Ref. [14J for more details). The non-local PNJL model can indeed be "derived" from 
QCD as demonstrated in Refs. p^6l [T7] . 

Once the input is fixed at zero temperature by well-known properties of the pseudoscalar mesons, 
the thermodynamics of the PNJL model can now be investigated with focus on the symmetry breaking 
pattern and on the intertwining of chiral dynamics with that of the Polyakov loop. The primary role 
of the Polyakov loop and its coupling to the quarks is to suppress the thermal distribution functions of 
color non-singlets, i.e. quarks and diquarks, as the transition temperature T c is approached from above. 
Color singlets, on the other hand, are left to survive below T c . This is seen by analyzing the relevant 
piece of the thermodynamic potential Q = —(T/V) \nZ ) the one involving the quark quasiparticles: 

AO ^ I ePp ^ r + ^ e _ (Bp _ M)/T + 3^ e -2(E p -»)/T + e -3(E p -»)/T 



T J (27T) 

d 3 p 

W) 



+ [ l n [i + 3 $* e -^)/ T + 3$ e ~ 2 ^ + ^ T + e -^ E " + ^ T ] . (5) 
J \2tiY l - 1 



This suppression of color non-singlets in the hadronic phase below T c should, however, not be interpreted 
as dynamical confinement. The unsuppressed color singlet three-quark degrees of freedom are not 
clustered but spread homogeneously over all space. At this stage, nucleons are not treated properly as 
localized, strongly correlated three-quark compounds plus sea quarks. Thus one cannot expect that such 
a model describes properly the low-temperature phase of matter at finite baryon chemical potentials 
around 1 GeV, the domain of nuclear many-body systems. 

At zero baryon chemical potential, a remarkable dynamical entanglement of the chiral and decon- 
finement transitions is nonetheless verified in the PNJL model, as demonstrated in Figj3} In the absence 
of the Polyakov loop the quark condensate (left dashed line) , shows the expected chiral crossover tran- 
sition, but at a temperature way below and far separated from the 1st order deconfinement transition 
controlled by the pure-gauge Polyakov loop effective potential U (right dashed line) . Once the coupling 
of the Polyakov loop field to the quark density is turned on, the two transitions move together and 
end up at a common transition temperature around 0.2 GeV. The deconfinement transition becomes a 
crossover (with Z{3) symmetry explicitly broken by the coupling to the quarks). 

It is interesting to compare this PNJL description of the interplay between the chiral and decon- 
finement transitions, with a scenario based on strong-coupling lattice QCD with inclusion of Polyakov 
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Figure 3: Solid curves: temperature dependence of the chiral condensate (qq) = (uu) + (dd) (left) and 
of the Polyakov loop <£> (right) normalized to the transition temperature T c = 205 MeV, as calculated in 
the nonlocal PNJL model with Nf = 3 and physical quark masses [15]. The dashed lines show the chiral 
condensate for the pure fermionic case and the Polyakov loop for the pure gluonic case, respectively. 



loop dynamics [19J. Although the starting points and frameworks are quite different in those two ap- 
proaches, recent studies [19] arrive at similar results concerning the close correspondence between chiral 
and deconfinement transition temperature ranges. 

Pions also have a pronounced influence [14] on the chiral condensate as it approaches the transition 
region from the hadronic phase. The chiral crossover becomes significantly smoother close to T c . Lattice 
QCD results [2j with physical pion masses and realistic extrapolations to the continuum limit show a 
similar tendency. 

Undoubtedly a prime challenge in the physics of strong interactions is the exploration of the QCD 
phase diagram at non-zero baryon density, extending from normal nuclear matter all the way up to 
very large quark chemical potentials ji q at which color superconducting phases are expected to occur. 
PNJL calculations at finite \± q give a pattern of the chiral order parameter in the (T, ji q ) plane showing 
a crossover at small \i q that ends at a critical point. From there on a first-order transition line extends 
down to a quark chemical potential \i q ~ 0.3 - 0.4 GeV at T = 0. 

The typical phase diagram that emerges from the most recent version of the non-local three-flavor 
PNJL model [15] is shown in FigJIJ The phase to the right of the 1st order transition line, but below the 
deconfinement boundary that decouples from the chiral transition beyond the critical point, has been 
named "quarkyonic" [21] and is hypothetically assumed to exist until superconducting phases take over 
at large quark chemical potential. A recent discussion about the possible relationship between such a 
hypothetical phase and the hadron rates produced in high-energy heavy-ion collisions [22] may not yet 
be conclusive but is pursued with great intensity. 

One must recall again at this point that calculations and extrapolations based on PNJL or related 
approaches still have limited validity at moderate and high baryon densities since they do not properly 
incorporate baryonic degrees of freedom. They miss the important constraints imposed by the existence 
of nucleons and nuclear matter. This becomes clear when converting the PNJL phase diagram, FigJIJ 
into a plot that translates quark chemical potentials into baryon densities via the derivative of the pres- 
sure, Pb = OP/ djiB ? with respect to baryon chemical potential \ib = 3/i g . The 1st order transition line 
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Figure 4: Phase diagram for the three-flavor nonlocal PNJL model [15j. The orange and blue bands show 
the deconfinement and chiral crossover transitions. The solid black line indicates the chiral first-order 
transition. 



converts to a broad coexistence region along the density axis that encloses what is well known as nuclear 
physics terrain, but now described with the "wrong" quasiparticle degrees of freedom: PNJL quarks 
rather than interacting nucleons. Quarks are not the relevant active quasiparticles at low temperature 
and baryon chemica potentials. Consequently, interpretations based on the schematic phase diagram of 
Figj4] along the chemical potential axis should be exercised with caution. 

Several further important questions are being raised in this context. A principal one concerns the 
existence and location of the critical point in the phase diagram. Extrapolations from lattice QCD, either 
by Taylor expansions [4j around \i = or by analytic continuation from imaginary chemical potential 
[23] . have so far not reached a consistent conclusion. A second question relates to the sensitivity of 
the first order transition line in the phase diagram with respect to the axial U(1)a anomaly in QCD 
[24] . This issue has been addressed in Ref. [25] using a general Ginzburg-Landau ansatz for the chiral 
SU(3) effective action with axial U(l) anomaly. It was pointed out that, depending on details of the 
axial U(1)a breaking interaction, a second critical point might appear such that the low-temperature 
evolution to high density is again just a smooth crossover, or the first order transition might disappear 
altogether and give way to a smooth crossover throughout. In a three-flavor NJL type realization of 
this approach [26], it has recently been demonstrated how the interaction between chiral condensate 
and diquark degrees of freedom, based on the genuine U(1)a breaking six-point vertex, can open up a 
smooth transition corridor along the chemical potential axis in which chiral and diquark condensates 
may coexist. 

Interesting developments are presently conducted towards fixing the input of the non-local PNJL 
model by direct comparison with lattice QCD thermodynamics at imaginary chemical potential [23], [271 
128] and then studying implications on the phase diagram at real \i. 

2.3 The Polyakov - quark-meson model 

A variant of the PNJL approach is the Polyakov - quark-meson (PQM) model [29j . It is based on a linear 
sigma model with quarks coupled to chiral pion and sigma fields and propagating in a Polyakov loop 



7 



200 



> 




100 



150 



50 




X first order 











100 



200 



300 



Mq[MeV] 



Figure 5: Phase diagram for the two-flavor PQM (Polyakov - quark-meson) model [30J. The solid and 
dotted lines show the deconfinement crossover transition (<£> and <3> crossover), while the dashed line 
marks the chiral crossover transition line (x crossover) ending in the critical point. The chiral first order 
transition is restricted to a small low-temperature region. 



background. Combined with renormalization group techniques, the PQM model has recently improved 
in a version [30j that includes the back-reaction of the quarks at finite \i q on the Polyakov loop effective 
potential. With increasing quark chemical potential, this back-reaction has the effect of aligning the 
chiral and deconfinement transition borders in the T — \± diagram such that these transitions tend to 
go in parallel not only at small chemical potential, but also as ji q increases. This important feature is 
apparent in Figj5} Note that the critical point is now moved to low temperature and the chiral first 
order transition is restricted to a very small region around quark chemical potential \i q ~ 0.3 GeV in 
the phase diagram. However, once again, this region corresponds to the density range of nuclear matter 
where the quark-meson picture is not expected to be valid. 

2.4 Dyson-Schwinger approach to QCD thermodynamics 

The approach that is, by design and construction, closest to QCD itself is the one based on solutions 
of (truncated) Dyson-Schwinger equations for quark and gluon correlation functions. Dyson-Schwinger 
thermodynamics has now advanced to the point of producing its generic version of the QCD phase 
diagram [31j, as displayed in Figj6| It is again remarkable that the chiral and deconfinement crossovers 
go jointly in parallel from T c ~ 170 — 180 MeV at ji q = down to a critical point at T e ~ 100 MeV and 
/jiq ~ 280 MeV. From there on downward a chiral first order transition region is indicated. 

What the Dyson-Schwinger and PQM approaches have in common, however, with the PNJL models 
is the absence of nucleon formation in the range of densities characteristic of nuclear (rather than 
constituent quark) matter. The range of quark chemical potentials around \± q ^ 0.3 GeV and low 
temperatures is nuclear terrain! The baryon densities 



determined by the derivative of the pressure P with respect to chemical potential, turn out to be in the 




(6) 
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Figure 6: Phase diagram derived from the Dyson- Schwinger approach to QCD thermodynamics [3T] . 
The curves show the chiral and deconfinement crossover transitions that turn out to almost coinicide 
over the whole phase diagram. The chiral crossover ends at a critical point located at T ~ 100 MeV 
and around \i q ~ 0.3 GeV, below which the chiral transition becomes first order. 



range 0.1 — 0.2 fm, covering typical nuclear matter densities. Therefore models and approaches dealing 
with dynamical (constituent) quarks as quasiparticle degrees of freedom in that area cannot be correct: 
it is crucial to take into account the clustering of quarks into nucleons and their correlations in the 
nuclear medium. 



3 Nuclear chiral thermodynamics 

At this point it is now appropriate to turn towards the discussion of nuclear thermodynamics using 
the framework of chiral effective field theory, the low-energy realization of QCD in its sector with 
spontaneously broken chiral symmetry. As a starter, consider the characteristic scales at work in the 
ground state of symmetric nuclear matter, i.e. with equal number of neutrons and protons (N — Z), 
illustrated in FigjTJ 

The nuclear many-body problem displays several "small scales" compared to the spontaneous chiral 
symmetry breaking scale, 4^/^ ~ 1 GeV. The Fermi momentum, p F ~ 1.36 fm -1 , is about half the in- 
verse pion Compton wavelength, A^ = m" 1 . The corresponding average distance between two nucleons, 
cInn — 1.8 fm, is about 1.3 A^. The nuclear compression modulus, k = (260 ± 30) MeV, is about twice 
the pion mass. All these scales are "pionic", i.e. of the same order as the pion Compton wave length 
or its inverse. The only exception is the binding energy per nucleon, |i?|/^4 ~ 16 MeV, an abnormally 
small number that reflects fine-tuning between short-range repulsion and intermediate-range attraction 
in the nucleon-nucleon force. 

With these observations in mind, it is suggestive that the long- and intermediate-range behavior of 
the effective interactions responsible for properties of nuclear matter around its equilibrium density can 
be described by a chiral expansion (in-medium chiral perturbation theory) in powers of pp <C ^nf^ ~ 1 
GeV. Short-distance correlations reflect large nucleon-nucleon scattering lengths and should instead be 
incorporated non-perturbatively by appropriate resummations or through a few contact interactions to 
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Figure 7: Illustration of the position of nuclear matter in the QCD phase diagram, outlining relevant 
scales of the nuclear many-body problem. 



be fitted to selected empirical data. 

Within the last decade such an approach to the nuclear many-body problem has been developed 
based on the understanding of the nucleon-nucleon interaction itself in terms of chiral effective field the- 
ory [32l ESI 133] ■ I n this approach, chiral one- and two-pion exchange processes in the nuclear medium 
are treated explicitly while unresolved short-distance dynamics are encoded in contact interactions. 
Three-body forces emerge naturally and play a significant role in this framework. The pion mass m^, 
the nuclear Fermi momentum p F and the mass splitting Ma — ~ 2m^ between the A (1232) and 
the nucleon are all comparable "small" scales that figure as expansion parameters. The relevant, ac- 
tive degrees of freedom at low energy are therefore pions, nucleons and A isobars. Two-pion exchange 
interactions produce intermediate-range Van der Waals - like forces involving the large spin-isospin po- 
larizablity of the individual nucleons. The Pauli principle plays an important role acting on intermediate 
nucleons as they propagate in two-pion exchange processes within the nuclear medium. 

3.1 In-medium chiral perturbation theory 

A key ingredient of chiral perturbation theory applied to nuclear matter is the in-medium nucleon 
propagator, 

Sn(p) = (^ + M n ) 

that takes into account effects of the filled nuclear Fermi sea (with Fermi momentum pp). These prop- 
agators enter in the loop diagrams generating the free energy density. Thermodynamics is introduced 
using the Matsubara formalism. The free energy density is then computed as a function of temperature 
T and baryon density p at given Fermi momentum p F . It is derived in the form reminiscent of a virial 
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(7) 
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expansion: 



roc roc roc 

T(p,T) = 4 dppKi(p)D{p) + / d Pl dp 2 K 2 (p 1 ,p 2 )D(p 1 )D(p 2 ) 
Jo Jo Jo 



+ 



dpi 



roc 

dp 2 / rfp 3 /C 3 (pi,p 2 ,P3) D(p 1 )D(p 2 )D(p s ) + ... (8) 
Jo 



Complete expressions for the kernels /Q are given in Refs.[34l 135] . For example, the one-body kernel 
/Ci represents the contribution of the non-interacting nucleon gas, including relativist ic corrections, to 
the free energy density and it reads 



/Ci M N + p 



pi 



3M t 



N 



4 
Pi 

8M 3 N ' 



(9) 



The terms involving K 2 and /C3 describe the effects of interactions and include two- and three-nucleon 
correlations, and so forth. The quantity 



1 + exp 



p 2 /2M N - fi 
T 



(10) 



denotes the density of nucleon states in momentum space, written as a product of the temperature 





Figure 8: One- and two-pion exchange processes 
contributing to the ground state energy (or the free 
energy density at non-zero temperature) . Each nu- 
cleon line is subject to a medium insertion (see 
Eq.@). 



Figure 9: Two-pion exchange processes involving 
one or two A (1230) isobars in intermediate states. 
The lower part of the figure shows the correspond- 
ing 3-body interaction induced by virtual A (1230) 
excitation. 



dependent Fermi-Dirac distribution and a kinematical prefactor p/27i 2 that has been included in D(p) 
for convenience. stands for the (free) nucleon mass. The baryon density p is calculated as 

roc 

p = 4, dppDip). (11) 
^0 

This relationship determines the dependence of the effective one-body "chemical potential" /i(p, T; Mn) 
on the thermodynamical variables (p, T). The "true" nucleon chemical potential that figures as a 
variable in the phase diagram, is determined by p = F + p(dF/dp) where F = T j p is the free energy 
per nucleon. The pressure is calculated using the standard relation 

P(ft7W ^. (12) 
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Figure 10: Equation of state of symmetric (TV = Z) nuclear matter from in-medium chiral effective field 
theory [34|,|35]. The critical point for the first-order liquid-gas phase transition appears at T c ~ 15 MeV. 
The equlibrium point at T = and baryon density p — 0.16 fm -3 is fixed by a single parameter, the 
strength of an NN contact interaction. All remaining input is pre-determined by known pion-nucleon 
interactions in vacuum. 



The loop expansion of in-medium chiral perturbation theory is related in a one-to-one correspondence 
to a systematic expansion of the energy density in powers of pp, with expansion coefficients expressed 
as functions of the dimensionless ratio vn^jpp. Calculations of nuclear thermodynamics have been 
performed up to and including three-loop order in the energy density. A typical set of pion-exchange 
diagrams is displayed in Figj8| Also explicitly included in the calculation are intermediate virtual 
A(1232) excitations and resulting three-body terms shown in Figj9j The three-loop truncation in 
the chiral expansion of the pressure in powers of the Fermi momentum, p F) valid unless four-nucleon 
correlations are substantial, implies that these calculations can be safely trusted up to about twice the 
density of normal nuclear matter, i.e. for p F less than 0.3 GeV <C 4^/^, with = 0.09 GeV the pion 
decay constant in vacuum. 

A limited set of constants associated with NN contact terms is fixed by reproducing e.g. the 
binding energy per nucleon in equilibrium nuclear matter at T = 0. Then this scenario leads to a 



realistic nuclear matter equation of state [34] (see Fig. 10) with a liquid-gas first order phase transition 
and a critical temperature of about 15 MeV, close to the range of empirical values extracted for this 
quantity. Incidentally, this nuclear liquid-gas transition is so far the only well established part of the 



phase diagram for strongly interacting matter at finite density. Fig. 11 shows the (T, ji B ) phase diagram 



for N = Z nuclear matter, while Fig, 12 displays (T, p) phase diagrams for isospin-asymmetric nuclear 
matter, outlining the evolution of the critical point and the liquid-gas coexistence region for increasing 
proton fractions Z/A [35]. The first order phase transition stops existing at Z/A = 0.05, at which point 
the liquid component disappears and neutron matter begins to be realized as an interacting Fermi gas. 
It is important to note that this behaviour is almost entirely driven by the isospin dependence of the 
in-medium two-pion exchange interactions. Typical ingredients of one-boson exchange phenomenology 
(such as rho and sigma bosons) do not appear in this approach and are replaced by explicit 2tt exchange 
mechanisms. 

Extensions of this framework to finite nuclei, using the energy density functional formalism, have 
been applied successfully [36] to calculate ground state properties (binding energies, mean square radii, 
deformations and systematics along isotopic chains) throughout the nuclear chart, from 16 O to Pb 
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Figure 11: Phase diagram of symmetric (TV = Z) 
nuclear matter from in-medium chiral effective 
field theory [35] . 



Figure 12: Evolution of the liquid-gas coexistence 
region from symmetric nuclear matter (Z/A = 0.5) 
to neutron-rich matter calculated using in-medium 
chiral effective field theory [35] . 



isotopes. Other examples of detailed nuclear structure phenomena such as the anomalously weak 
14 C beta decay transition, with its lifetime of almost six thousand years that enables radiocarbon 
dating, can be understood in terms of the density dependent chiral effective interaction including a 
pronounced three-body contribution [37J . Further recent applications include the investigation of chiral 
NN interactions in the context of nuclear energy density functionals and the quasiparticle interaction 
in Fermi liquid theory [38] . 



3.2 Density and temperature dependence of the chiral condensate 

In a nuclear equation of state based on chiral dynamics the pion mass enters explicitly (or, equivalently, 
the quark mass m q according to the Gell-Mann - Oakes - Renner relation, m^fl = -m g (^)). Equiped 
with such an equation of state one can now ask the following questions: how does the chiral (quark) 
condensate extrapolate to baryon densities p exceeding those of normal nuclear matter, and how 
does it evolve as a function of temperature T in a nuclear environment? The answer, first at T = 0, 
is based on the Hellmann-Feynman theorem applied to the Hamiltonian density of QCD in the nuclear 
ground state |^), 



31-1 

WVI*> = <*|-^^l*> = <W( 



fn dm l 



(13) 



where (^)o is the chiral condensate in vacuum. The energy density £ (normalized so that it vanishes 
at p — >> 0) is expressed at given baryon density p as a function of the pion mass. In-medium chiral 
effective field theory gives the following result [39j : 

W)p = i _ P a N L 3p| \ p d ( E int {m^p F ) \ (u . 

(W>)o flml\ 10M* N ^'") ftdml\ A J [ } 

The second term on the r.h.s., with its leading linear dependence on density, is the contribution from 
a free Fermi gas of nucleons. It involves the pion-nucleon sigma term a N ~ 0.05 GeV and additional 
non-static corrections. The third term includes the pion mass dependence of the interaction energy 
per nucleon, E int /A. This term features prominently the two-pion exchange interaction in the nuclear 
medium including Pauli principle corrections, and also three-nucleon forces based on two-pion exchange. 
The dashed-dotted curve in Fig fl~3] shows the pronounced leading linear reduction in the magnitude of the 
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Figure 13: Density dependence of the chiral condensate at temperature T = in symmetric nuclear 
matter [39]. Dot-dashed curve: leading order term using a N = 50 MeV. Upper curve: full in-medium 
chiral effective field theory result at three-loop order. Lower curve: chiral limit with vanishing pion 
mass. 



chiral condensate with increasing density. This holds in the absence of correlations between the nucleons. 
Up to about the density of normal nuclear matter, this term dominates, whereas the interaction part 
tends to delay the tendency towards chiral restoration when the baryon density is further increased. 
This stabilizing effect arises from the combination of Pauli blocking in two-pion exchange processes 
and three-body correlations, and their explicit dependence on the pion mass. This behaviour is highly 
sensitive to the actual value of the pion mass. In the chiral limit, — >> 0, with stronger attraction in 
the NN force at intermediate ranges, the trend is reversed and the more rapidly dropping condensate 
would now lead to the restoration of chiral symmetry at relatively low density. In essence, nuclear 
physics would look radically different if pions were exactly massless Nambu-Goldstone bosons. The 
drastic influence of explicit chiral symmetry breaking in QCD through the small but non-zero u and d 
quark masses on qualitative properties of nuclear matter is quite remarkable. 

The additional temperature dependence of the in-medium chiral condensate [40] is found by repeating 



the calculation, Eq.(13), using the free energy density, T{m^\ p,T): 



(15) 



The result shown in Fig|T4] demonstrates how the chiral condensate "melts" with increasing temperature 
and approaches the linear density dependence of the non-interacting Fermi gas at a temperature of about 
100 MeV. At low temperatures, the influence of the nuclear liquid-gas first-order phase transition is 
visible even in the chiral condensate. This becomes particularly pronounced when the condensate is 



plotted as a function of baryon chemical potential fi (see Fig. 15). The effect of thermal pions on the 
temperature dependence of the quark condensate must still be added; it amounts [HI [42] to lowering 
the magnitude of the condensate by another 5% of its vacuum value at 100 MeV. 

However, no tendency towards a first order chiral phase transition is visible, at least up to baryon 
densities about twice that of normal nuclear matter and up to temperatures T ~ 100 MeV. In that 
range it appears that the relevant fermionic quasiparticle degrees of freedom are indeed nucleons rather 
than quarks. This should be kept in mind as an important constraint for the discussion of the QCD 
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Figure 14: Density and temperature depen- 
dence the chiral condensate in symmetric nu- 
clear matter [35]. At T = the influence of 
the first-order liquid-gas phase transition is in- 
dicated by drawing a straight line corresponding 



the Maxwell construction (ignored in Fig, 13) 



Figure 15: Temperature dependence of the chi- 
ral condensate [35j in symmetric nuclear matter 
as function of the baryon chemical potential \i. 
This plot demonstrates the pronounced influ- 
ence of the first-order liquid-gas phase transi- 
tion at low temperatures. 



phase diagram at finite baryon chemical potential. 



4 New constraints from neutron stars 

The equation-of-state of highly compressed cold matter in the center of neutron stars has been an issue 
and a challenge for many decades. For a long time it was thought that the vast majority of neutron 
stars have masses around 1.5 times the solar mass. The equation-of-state of neutron star matter that 
serves as an input for the Tolman-Oppenheimer-Volkov equation to calculate the mass-radius relation 
of the star, M{R) ) should then be sufficiently "stiff" in order to support such neutron star masses. 
Many scenarios for building the cores of such compact stars were designed and discussed (43] 5 from 
conventional nucleonic matter equations-of-state [44] to matter with more exotic components such as 
kaon condensates and quark matter [45j . 

The empirical constraints were considerably sharpened recently when a neutron star in a binary with 
a white dwarf was observed. An extremely accurate Shapiro delay measurement of the pulsar signal 
resulted in the unusually large mass of 1.97 ± 0.04 in solar mass units [46]. This observation rules out 
some of the softer equations-of-state, e.g. those featuring kaon condensation. 

A further important constraint comes from the systematic analysis of a well examined series of 
neutron stars. Using this data base, a most probable mass-radius relation M(R) has been established 
[47] . with a relatively narrow window of star radii around R ~ 12 km. The combination of all these 
pieces of information can now be used in order to set rather tight constraints for the equation-of-state of 
dense neutron star matter. One starts out at low densities with a reliable equation-of-state for the crust 
of the star and then proceeds to higher densities using a realistic neutron matter equation-of-state based 
on the chiral effective field theory calculations discussed earlier [34j [35J [48] . Such an extrapolation holds 
up to baryon densities of about twice that of normal nuclear matter. The further extrapolation to high 
densities is then performed using sections of polytropes, parametrizing the pressure as P = const. p r . 
When this is all done one arrives at an equation-of-state, i.e. pressure P = P(S) as a function of energy 
density £, with an uncertainty band determined by the new empirical constraints [48l [49] . The result of 
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Figure 16: Pressure versus energy density for neutron star matter (uncertainty band included) [49] 
with built-in constraints resulting from the existence of a two-solar-mass neutron star [46] and from 
the analysis of most probable mass-radius relations, M(i?), from Ref. [47]. The black solid area is the 
neutron matter equation of state (with 10% protons) calculated using in-medium chiral effective field 
theory to three-loop order [35]. The solid line is the APR equation of state [44] based on a realistic 
nuclear many-body calculation. 



such an analysis is shown in Fig(l6} The condition of supporting a two-solar-mass neutron star implies 
that this equation-of-state must be quite stiff, with a high-density adiabatic index T~ 3.7 — 3.8. It is 
remarkable that, at present, the only equation-of-state that meets all constraints mentioned is the one 
based on standard nucleon degrees of freedom [44] with realistic nuclear two- and three-body forces. 
For a more detailed discussion of the implications of such analysis, see Ref. [50] . 



5 Summary and Outlook 

From the variety of existing model calculations (including those using PNJL models) one might draw 
the presumably premature conclusion that critical phenomena in compressed baryonic matter occur 
already at a density scale not much higher than that of normal nuclear matter. However, these models 
are so far not capable of working with the proper degrees of freedom around and below baryon chemical 
potentials ~ 1 GeV (corresponding to quark chemical potentials around 0.3 GeV). Approaching the 
corresponding density scale from below, it is obvious that constraints from what we know about the 
nuclear matter equation of state must be seriously considered. The implementation of such nuclear 
matter constraints (that are well reproduced using the framework of in-medium chiral effective field 
theory) into a consistent picture that combines the hadronic sector of the QCD phase diagram with the 
high-density and high-temperature regions remains a major challenge for the future. 

Finally, in the discussion of high density, low temperature constraints on the QCD equation of state, 
recent new developments concerning the mass-radius relation deduced from observations of neutron 
stars in binaries are of great interest in reducing the set of acceptable equations of state. In particular, 
too soft exotic equations of state are most likely ruled out by the existence of a two-solar-mass neutron 
star, whereas an equation of state with "standard" nuclear physics degrees of freedom and realistic 
correlations including three-body forces appears to be consistent with observations. 
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